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Shuttle  Contamination  Modeling:  The  Plasma 
Wave  Field  of  Spacecraft 


I.  IIMTRODI'CTION 

Shuttle  and  other  low  altitude  Earth  orbit  spacecraft  move  through  a  medium 

9  *3 

that  is  characterized  by  a  neutral  number  density  of  about  10  cm  and  an  ion  den- 
4  5  -3 

sity  of  10  -  10  cm  .  Typical  collisional  mean-free  paths  are  on  the  order  of 

kilometers  or  greater,  precluding  a  strong  collisional  interaction  with  the  space¬ 
craft  which  would  produce  collision-dominated  waves  and  shocks.  The  long  mean- 
free  paths  also  mean  that  the  interaction  between  the  ambient  ions  and  neutrals  is 
weak  on  the  scale  of  kilometers,  so  that  the  ions  may  be  treated  as  a  collisionless 
plasma  on  this  scale.  Typical  temperatures  at  shuttle  altitude  are  on  the  order  of 
1000  K,  so  that  the  electron  thermal  speed  is  on  the  order  of  100  km /sec,  about 
an  order  of  magnitude  less  than  the  Alfven  speed.  As  a  result,  the  plasma  may  be 
treated  in  the  cold  plasma  approximation.  A  spacecraft  moving  through  such  a 
cold,  collisionless  plasma  can  be  expected  to  produce  waves  and  shocks  or  shock 
like  structures. 

Typical  phase  and  group  velocities  of  cold  plasma  waves  are  on  the  order  of 
the  Alfven  speed.  As  such,  shock  waves  could  be  regarded  as  unlikely,  because 
the  spacecraft  velocity  is  a  few  kilometers  per  second,  quite  small  compared  to 
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the  wave  speeds.  However,  cold  plasma  waves  range  from  almost  isotropic  to 

quite  anisotropic.  For  some  wave  modes,  the  phase  or  group  velocity  goes  to 

1  2 

zero  in  certain  directions  with  respect  to  the  background  geomagnetic  field.  ’ 

The  possibility  exists  that  these  modes  can  be  shocked  by  any  spacecraft  velocity, 
no  matter  how  small. 

The  interaction  between  spacecraft  and  the  ambient  plasma  has  been  investi- 
3  4 

gated  by  Drell  et  al  and  by  Barnett  and  Olbe rt.  Drell  et  al  showed  that  space¬ 
craft  could  produce  Alfven  wings,  a  steady  Alfven  disturbance  extending  away  from 
the  spacecraft.  Barnett  and  Olbert  have  recently  extended  the  analysis  to  include 
all  plasma  waves,  although  a  written  report  is  not  yet  available. 

The  purpose  of  the  present  paper  is  to  describe  the  wave  and  shock  structure 
of  a  spacecraft  moving  through  a  collisionless  cold  plasma.  The  focus  of  the  work 
will  be  the  spatial  properties  of  those  waves,  produced  by  the  interaction  of  the 
spacecraft  with  the  ambient  plasma,  which  can  propagate  to  great  distances  from 
the  spacecraft  and  can  form  shock-like  disturbances.  These  waves,  unlike  waves 
which  propagate  isotropically,  can  propagate  only  within  sharply  defined  spatial 
regions.  Our  methods  lead  to  definite  predictions  for  (1)  the  frequencies  of  waves 
which  can  propagate  to  large  distances  from  the  spacecraft,  (2)  the  directions 
within  which  they  can  propagate,  (3)  the  wave  power  spectrum  around  the  space¬ 
craft,  and  (4)  predictions  of  shock-like  disturbances. 

An  important  restriction  of  our  results  should  be  mentioned  here.  We  use 
essentially  two-dimensional  methods  throughout.  The  formalism  involved  in  the 
three-dimensional  problem  is  sufficiently  complicated  that,  without  the  two-dimen¬ 
sional  formulation  as  introduction,  the  results  are  essentially  unreadable.  Quali¬ 
tatively,  a  two-dimensional  treatment  is  no  great  impediment.  The  typical  proper¬ 
ties.  exemplified  by  the  characteristics,  are  similar  in  three  dimensions.  How¬ 
ever,  the  quantitative  details  should  be  regarded  and  used  with  a  certain  amount  of 
discretion;  the  details  can  be  expected  to  change  in  three  dimensions. 

One  motive  for  this  work  is  to  begin  a  fundamental  investigation  of  the  inter¬ 
action  of  spacecraft  with  the  plasma  environment.  Such  an  investigation  is  neces¬ 
sary  to  treat  a  wide  range  of  issues,  some  of  traditional  importance,  some  newly 
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emerging.  These  issues  can  he  broadly  classed  under  two  main  headings:  space¬ 
craft  charging  and  spacecraft  contamination. 

Spacecraft  charging  issues  include  charging  of  the  spacecraft  itself,  wake 
formation,  and  density  enhancements  in  ram.  These  are  all  affected,  at  least  in 
principle,  by  plasma  waves  and  turbulence.  A  particularly  pressing  issue  is  the 
interaction  of  the  wake  with  waves  and  turbulence  and  the  consequent  effect  on 
charging. 

Spacecraft  contamination  issues  include  spacecraft  glow,  ionization  of  high 
density  contaminants  of  spacecraft  origin,  and,  possibly,  extended  plasma  radia¬ 
tion  processes.  While  the  present  work  cannot  hope  to  address  all  of  these  topics, 
it  should  at  least  provide  the  basis  for  addressing  several  of  them,  including  ef¬ 
fects  of  waves  on  charging  and  the  evolution  of  plasma  clouds. 

A  second  motive  for  this  work  is  to  provide  the  mathematical  and  physical  ba¬ 
sis  for  modeling  of  the  spacecraft  interaction.  Current  efforts  are  largely  limited 
to  purely  numerical  work  (for  example,  the  NASCAP  and  POLAR  codes!  because 
of  the  intractability  of  the  physics  as  well  as  the  complications  arising  from  the 
geometry  of  real  spacecraft.  These  codes  are  typically  unwieldy  and  limited  by 
the  simplifying  assumptions  necessary  to  produce  a  computer  program  that  runs  in 
a  finite  length  of  time.  If  one  contemplates,  as  we  do,  computer  modeling  of  a 
wide  range  of  spacecraft  interactions,  including  charging,  contamination,  surface 
chemistry,  chemical  reactions,  molecular  transport,  ionization  processes,  and 
the  formation  of  plasma  waves  and  turbulence,  while  perhaps  including  the  effects 
of  charged  particle  beams,  the  probability  of  success  with  a  purely  numerical  ap¬ 
proach  appears  to  be  negligible:  the  reason  one  cannot  see  the  light  at  the  end  of 
the  tunnel  is  because  there  is  no  end  to  the  tunnel.  The  approach  that  we  advocate, 
and  that  has  had  reasonable  success  in  the  development  of  the  NASCAP  and  POLAR 
codes,  is  to  reduce  each  part  of  the  problem  to  manageable  proportions  using  ana¬ 
lytical  techniques  where  possible.  This  paper,  we  hope,  represents  a  contribution 
to  that  process. 

In  the  following  sections,  we  will  treat  cold  plasma  waves  in  a  homogeneous 
medium  directly  from  the  point  of  view  of  the  governing  differential  equations. 
From  the  equations,  we  will  determine  the  characteristics  that  govern  wave  prop¬ 
agation  and  show  that  there  are  three  different  frequency  domains  in  which  waves 
can  propagate  to  large  distances  from  the  spacecraft.  The  propagation  of  discon¬ 
tinuities  is  discussed  in  terms  of  and  related  to  the  characteristics.  The  relation 
between  cold  plasma  resonances  and  the  characteristics  of  the  differential  equa¬ 
tions  is  pointed  out,  and  a  physical  picture  of  the  characteristics  as  carriers  of 
discontinuities  is  developed  from  the  plane  wave  picture.  Solutions  of  the  cold 
plasma  equations  appropriate  for  localized  disturbances  are  then  determined  ana- 


lytically;  they  provide  a  prediction  of  the  power  spectrum  of  fluctuations  far  from 
the  spacecraft.  Finally,  applications  to  the  spacecraft  environment  are  discussed. 


2.  COLD  PLASMA  WAVES 

The  treatment  of  plasma  waves  begins  with  a  standard  approach  to  the  prob¬ 
lem.  We  include  it  to  establish  approach  and  notation.  The  aim  is  to  write  partial 
differential  equations  governing  waves  in  a  cold  plasma.  In  particular,  we  shall 
reduce  the  problem  to  a  fourth  order  partial  differential  equation  for  a  magnetic 
field  fluctuation. 

The  linearized  cold  plasma  equation  of  motion  is 

(d!dt)6\a  =  (e^/m^ldE  +  (ea/mac)5Va  x  B  (1) 

where  «  indicates  the  species  of  particle,  ea  is  its  charge,  ma  its  mass,  dVa  is 
the  bulk  velocity  of  the  species,  ^E  is  the  electric  field,  6  indicates  fluctuating 
quantities,  and  c  is  the  speed  of  light.  The  solution  for  iV  with  a  harmonic  time 
dependence,  did t  =  -icu,  is 

<5 V„  =  (ealma){iulu2)  <5E  (2) 

4^  =  (eal  ma)(iu>6Ex  +  «5E  x  Qa)  /  (J  -  Q2)  (3) 

where  ||  and  ^  mean  parallel  and  perpendicular  to  the  background  magnetic  field, 
and  the  gyrofrequency  of  species  a  is  defined  by 

Qa  =  -eaB/ (mac).  (4) 

We  note  that  the  frequency  W  is  defined  in  the  proper  frame  of  the  plasma  (the 
frame  in  which  the  background  plasma  is  at  rest)  and  that  the  frequency  measured 
by  an  observer  moving  through  the  plasma  with  a  velocity  U  is  u'  =  u  *  k*U,  where 
k  is  the  wave  vector. 

The  current  density,  defined  by 

6  J  -■  EanaeadV  (5) 

where  na  is  the  number  density  of  species  a,  is 


iJ  =  (iu/47T)(wp2/tJ2)  dE|j  +  (icj/47T)i;a[Up^2/((J2  -  .Qa2n^Ei- 


+  <  1  / 4 rr ) [Z^Wp cx2i20t / <^2  -  Oa2)}b  x  <5  E 


(6) 


where 


upa2  =  4rrnaea2/ma 

(7) 

is  the  plasma  frequency  of  species  a,  the  plasma  frequency  is 

2  v'  2 

oj  -Y  cj 

p  Ot  pOt 

(8) 

A  ♦» 

and  b  =  B/B  is  a  unit  vector  in  the  direction  of  the  field.  Inserting  Eq. 

(6)  into  the 

combined  Ampere -Faraday  Law 

P*<?  x  dE)  =  (47rio»/c2)i"j  +  (u2/c2)iE 

(9) 

we  obtain 

PxiPxiE)--  f0iEll  "  iE 

(10) 

where 

(o  =  (w2  -  V)/c2 

(11) 

€P=  -{(J2lc2)  [1  -Zupa2l(u,2  -  Q2)] 

(12) 

«H  =  -<iw/c2)i:Wpc,2fia/(w2  -  fi«2)- 

(13) 

Our  approach  will  be  to  write  the  differential  equations  in  a  homogeneous 
background  plasma.  To  do  this,  we  assume  that  the  magnetic  field  is 


H  -  (0,  C,  B  ) 


(14) 


a/ay  s  o 


(15) 


and  that  the  coefficients  €  are  constant.  Eq.  (15)  limits  the  subsequent  analysis  to 
two  dimensions;  this  approach  sacrifices  generality  for  the  sake  of  clarity.  The 
assumption  of  constant  magnetic  field  is  consistent  with  the  nature  of  the  space¬ 
craft  interaction;  the  spacecraft  is  a  small  object  in  a  large  scale  (on  the  order  of 
one  Earth  radius)  magnetic  field.  The  assumption  of  constant  f ,  which  implies 
that  the  number  density  is  constant,  is  not  very  good;  the  density  changes  by  or¬ 
ders  of  magnitude  in  a  large  region  occupied  by  the  spacecraft  wake  and  perhaps 
in  a  smaller  region  in  ram.  As  a  result,  the  detailed  results  of  wave  propagation 
will  be  limited  to  a  region  outside  the  strong  density  gradients  near  'he  spacecraft. 
It  is  not  difficult  to  generalize  our  differential  equations  to  variable  density,  but 
the  possibility  of  solving  the  equations  analytically  appears  to  be  remote.  With 
these  assumptions,  Eq.  (10)  becomes 


-djd  <5e  -  d  6 E  )  =  -tn  <5e  -  <„  0E 
z  z  x  x  z  P  x  H  y 


-  id  2  *  3  2)  jr:  .  =  -  f  <5E  -  f  D  (5E 
x  z  y  H  x  P  y 


d  ( d  6  e  -  d  6  E  )  =  ( n  6  K 
X  Z  X  X  7  0  z 


where  d.  means  <5  c?x  .  (.'sing  the  Faradav  relation 

l  t 


a,«5nv  =  -iu4Bv  =  -c<a7(5Ex  -  dx6 k2> 


where  0  mean.-  d  d',  Kqs.  (16)  through  (18)  can  be  written 


( i OJ  '  e  <d?  ‘  fpdE  " 


«V  •  V  -  V*Kv  ;  eHdi:x 


(iw  <  i  d  <JP  f  b  F  . 

XV  it  7 


Kqs.  (20>  and  <22>  van  he  combined  to  yield 


[0/  -  f0’as'  -  q.]^Bv  1  €Haz5i:v 
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or,  because 


<e/iw)dzdEy  =  -tBx. 

[az2-(ep/€0)ax2-  €p]dBy=  £HdBx.  (25 

Similarly.  Eq.  (21)  carrbe  written  (after  differentiation  with  respect  to  z) 

<dx2  +  dz2-  «p>*Bx=-eH(l+  (o'ld*)iBy  <2« 

where  Faraday's  Law  has  been  used  to  eliminate  derivatives  of  <}E. 

For  future  reference,  we  should  like  to  state  that  Eqs.  (25)  and  (26)  can  be 
rewritten  in  the  form 

fa2  -  (tp/fo^x2  -  «P]«  =  <H<1  +  Vlax2>£  (2‘ 

(ax2+az2-  (P>0=-(Ha  <2t 

where  at  and  are  defined  by 

a  =  (1  +  €p/€0)  dBv  -  (€H/€0)<5Bx  ( 2 £ 

P  -  (1  f  <p/f0)<JBx  +  UH/€0)dBy.  (3C 

Eqs.  (25)  and  (26)  or  Eqs.  (27)  and  (28)  are  a  pair  of  coupled  second  order 
partial  differential  equations  that  govern  linear  waves  in  cold  quasineutral  plas¬ 
mas.  It  is  instructive  to  consider  the  low  frequency  limit  of  these  equations.  In 
the  limit  that  0,  it  is  easy  to  see  that 

fp  =  -(o»2/c2)(l  +  c2/Va2) 


where 


is  the  nonrelativistic  Alfven  speed,  and 


-  (47riw/cH)Eanaea  =  0  (33) 

because  of  the  quasineutrality  of  the  background  plasma,  which  we  here  assume. 

In  this  case,  Kqs.  (25)  and  (26)  uncouple,  and  can  be  written 


[' dy 2  -  Up l(0)dx2  '  =  0  (34) 

<ax2  az2  -  <P>  <5BX  =  0.  (35) 

These  can  tic  recognized  as  the  Alfven  mode  and  fast  mode,  respectively.  The 
Alfven  mode  propagates  approximately  along  the  field,  while  the  fast  mode  propa¬ 
gates  isotropically. 

Ilqs.  (25)  and  (26)  can  lie  reduced  to  a  single  equation: 

•  dz2- V 

+  £I|2  (1  f  fo’laX2)  *Ky  =  °  (36) 

Kq.  (36)  is  a  single  fourth  order  equation  governing  cold  plasma  waves. 


3.  I  I IK  (  II  \K  VC  I  KKISTICS  (KIVHRMMI  WAYK  PROPAGATION 

By  a  standard  mathematical  procedure,  the  characteristics  of  the  differential 
Kq.  (36)  can  be  determined  by  solving  a  first  order  partial  differential  equation 

Q  (V<i»  0  (37) 

(cf.  C’ourant  and  Hilbert,  632® ). 
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In  the  case  of  Eq.  (36),  the  characteristic  equation  is 


<*x2  +  ^z2>[^2  *  WO  =  °-  (38) 

There  are  two  solutions  for  characteristic  surfaces: 

<PX2  +  <t>2  =  0  (39) 

and 

+x  =  <(0/€P^z2-  <4°> 

Eq,  (39)  yields  only  imaginary  surfaces.  In  the  terminology  of  partial  differ¬ 
ential  equations,  it  means  that  the  problem  is  elliptic.  Physically,  it  means  that 
the  waves  associated  with  these  characteristics  propagate  more  or  less  isotropic¬ 
ally.  An  example  of  completely  isotropic  propagation  is  the  fast  mode  wave,  de¬ 
scribed  by  Eq.  (35).  Because  of  the  isotropic  propagation,  the  wave  amplitudes 

can  be  expected  to  fall  off  rapidly  from  a  confined  source.  The  falloff  for  fast 
2 

mode  waves  is  l/r  .  While  the  falloff  will  be  less  rapid  for  more  anisotropic 
wave  modes,  it  can  be  expected  to  be  significant. 

Eq.  (40).  on  the  other  hand,  has  real  or  imaginary  characteristic  surfaces 
depending  on  the  wave  frequency.  For  frequencies  for  which  the  characteristics 
are  real,  the  problem  is  said  to  be  hyperbolic.  The  waves  associated  with  real 
characteristics  do  not  propagate  isotropically;  they  are  subject  to  shocks  and  dis¬ 
continuities.  In  this  linear  problem,  the  shocks  or  discontinuities  propagate  only 
along  the  characteristics,  as  do  any  discontinuities  in  the  derivatives. 

Even  when  Eq.  (40)  has  real  characteristics,  there  are  real  and  imaginary 
characteristics  in  the  same  equation.  The  effect  of  this  coupling  will  be  examined 
later  when  we  obtain  explicit  solutions.  Except  to  the  extent  that  the  waves  asso¬ 
ciated  with  imaginary  characteristics  are  driven  by  waves  as  ociated  with  real 
characteristics,  they  will  not  be  treated  in  detail  in  this  paper.  Such  waves  would 
be  worth  a  separate  investigation. 

The  solution  of  Eq.  (40),  obtained  by  the  method  of  characteristics,  is  <p  = 
constant  along  the  surfaces  defined  by 

dx  =  ±  y (fpleQ)dz.  (41) 

The  coefficients  and  ( p  can  each  be  positive  and  negative  as  a  function  of  fre- 
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3 


Based  on  these  results,  the  frequencies  for  which  there  are  real  characteris¬ 
tics  can  be  determined  (Table  1).  We  have  tacitly  assumed  the  ordering 
2  2 

Q  <  <  u  ,  typical  of  the  low  Earth-orbit  environment.  There  are  three  re- 
e  p 

gions  that  can  support  shock-like  disturbances:  (1)  Alfven--below  the  ion  gyrofre- 
quency,  (2)  whistler- -between  the  lower  cutoff  frequency  and  the  electron  gyrofre- 
quency,  and  (3)  upper  hybrid --between  the  plasma  frequency  and  the  upper  cutoff 
frequency.  In  these  regions,  the  angle  between  the  characteristics  and  the  magne¬ 
tic  field  is 

B  =  ±tan  V(fp/f0)*  (44) 


The  angle  0  is  shown  in  Figure  1  for  a  plasma  typical  of  low-altitude  Earth 


Figure  1.  The  Angle  Between  the  Characteristics  and  the  Magnetic 
Field,  Given  as  a  Function  of  Wave  I  renuenev,  for  a  Hydrogen  Plasma 
With  n  =  lO’i'.'n"'1  and  H  0.  3G 


orbit:  a  hydrogen  plasma  with  n  =  10  cm  and  B  =  0.  3G.  In  order  of  increasing 
frequency,  the  three  curves  represent  the  characteristics  associated  with  the  Alf- 
ven,  whistler,  and  upper  hybrid  regions.  At  very  low  frequencies,  the  angle  is 
essentially  zero,  the  well  known  result  for  Alfven  waves:  Alfven  waves  propagate 
only  along  the  field.  The  angle  broadens  as  the  frequency  increases  and  becomes 
90°  at  the  ion  gyrofrequency.  This  behavior  repeats  above  the  angle  is  zero 

at  anc*  *ncreases  to  90°  at  the  electron  gyrofrequency.  Above  the  plasma  fre¬ 

quency,  the  behavior  is  reversed:  the  angle  is  90°  and  decreases  to  zero  ati?^^. 
The  discussion  of  the  behavior  of  the  waves  with  respect  to  these  angles  is  deferred 
to  a  later  section. 

I  Vrofutfiittion  of  Di'Coiili/iiiitic- 

The  above  discussion  of  the  existence  of  discontinuities  originated  with  the 
formal  mathematical  treatment  of  the  characteristics,  which  is  quite  abstract  at 
best.  The  present  section  is  designed  to  show  the  relation  of  the  characteristics 
to  the  discontinuities.  The  treatment  is  not  general,  but  is  much  more  direct  and 
physical. 

To  treat  discontinuities,  we  start  with  Eqs.  (16)  through  (18)  for  the  electric 
field  components.  It  follows  from  Eqs.  (16)  and  (18)  that 

3x(-‘paExf  VV+aZ(€o*Ez)  =  0-  (45) 

This  is  a  divergence  equation.  Discontinuities  can  be  treated  by  integrating  over 
a  Gaussian  pillbox  surface  enclosing  the  discontinuity  and  using  the  divergence 
theorem 

fV’  Ad'*x  ^n«Ada.  (46) 

For  a  jump  discontinuity,  Eq.  (46)  shows  that  the  normal  component  of  the  vector 
A  is  continuous: 

[A  ]  0.  (47) 

L  n 


The  bracket  notation  is  defined  by 

[Q  ]  s  Q2  "  Qj 


(48) 


where  Q  is  any  quantity  and  1  and  2  refer  to  the  two  different  sides  of  the  discon¬ 
tinuity. 


Applying  this  method  to  Eq.  (45),  we  obtain 
nxt’€Pdl'x  f  £H<5Ky3  <  n?[(0iE7J  =  0. 


(49) 


Assuming  that  the  coefficients  f  are  continuous,  we  have 

-"xV  '  nx£H  tiI:yl  f  n7.eotiKz^  =  °‘  <50) 

In  a  similar  manner,  it  follows  from  Faraday's  l.aw  (see  Eqs.  (19)  and  (24)  that 


n/[iKx3  •  nx  =  0 

(51) 

[i.;y]  =  °. 

(52) 

The  only  solutions  of  Kqs.  (50).  (51),  and  (52)  for  which  the  electric  field  discon¬ 
tinuities  do  not  vanish  are  given  by 


V£e  -  "/‘o  5  0 


"x  *^<VW 


(53) 


(54) 


The  normal  vector  is  perpendicular  to  the  surface  of  discontinuity,  and,  by 
i;<j.  (44),  perpendicular  to  a  characteristic  surface.  It  follows  that  any  discontin¬ 
uities  in  the  electric  field  must  lie  on  a  characteristic  surface.  It  further  follows 
that  there  are  no  discontinuities  at  all  unless  the  characteristics  are  real. 


rv.v 


A  A 
.•■A 


VAA. 


V-'.V.O. 
*  .  «  .1  .  « 

-.V-.v 


.3.2  anil  Eronp  \  clocilio 

Since  the  methods  developed  here  use  the  same  assumptions  as  a  more  ortho¬ 
dox  plane  wave  treatment  of  plasma  waves,  it  is  to  be  expected  that  the  important 
features  should  be  the  same.  This  is  the  case  in  several  respects.  In  particular, 
the  characteristics  can  be  identified  as  the  resonance  lines  of  the  plane  wave  pic¬ 
ture.  The  purpose  of  this  section  is  to  provide  a  physical  picture  of  the  charac¬ 
teristics  as  resonances  to  aid  in  the  development  of  the  idea  of  the  characteristics 
as  carriers  of  discontinuities. 

The  dispersion  relation  of  cold  plasma  waves  car.  be  obtained  from  any  of  the 
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principal  equations:  Eqs.  (16),  (17),  and  (18),  Eqs.  (27)  and  (28),  or  Eq.  (36).  It 
is  convenient  to  write  it  in  the  form 


D(tt.  k)  =  a/k4  +  b/k2  +  c  =  0  (55 

where  the  coefficients  are  defined  by 

a=  V6P2  +  ‘h2)  (56 

b  =  2c0€pcos2^+  (<0fp  *  fp2  -  «H2)sin2^  (57 

c  =  (qCOs^iI/  -  fpSin2^  (58 

and  is  the  angle  between  k  and  B. 


To  investigate  the  relation  of  resonances  to  the  characteristics,  we  inquire 
about  the  behavior  of  the  phase  and  group  velocities  when  the  wave  vector  is  per¬ 
pendicular  to  the  characteristics.  If  k  is  perpendicular  to  one  of  the  characteris¬ 
tics,  defined  by  Eq.  (44),  then 

tan2 lj/  =  1/tan2 6  =  (Ql  €p  (59 

so  that 

c  -  0  (5® 

which  shows  that  the  characteristics  are  resonance  lines  (see  Stix,  p.  14)  and 


b  "  €0(*0€P 


eH  )/(f0  +  fP1, 


Then  the  roots  of  the  dispersion  relation  are 


We  conclude  that  when  k  is  perpendicular  to  a  characteristic,  the  associated  wave 
vector,  given  by  Eq.  (62),  is  infinite  and  the  phase  velocity  is 

V  .  =  <J/k  =  0.  (64) 

Ph 

Let  us  pursue  the  same  line  of  thought  with  the  group  velocity.  The  group 
velocity  is  defined  to  be 

V„  =  -(dDldhndDldu).  (65) 

g 

The  numerator  is 

dD/dk  =  ek<aD/dk)  +  (t^/kHdD/dtif)  (66) 

where  the  unit  vectors  e^  and  e  are  parallel  and  perpendicular  to  k,  respectively 
Evaluated  at  c(ijr)  =  0,  so  that  k  is  perpendicular  to  a  characteristic. 


Vgk  :  2b/(k3dc/dw)  =  0 

(67) 

Vgtfr  =  -{dc/d^)l{kdcldu)  =  0 

(68) 

so  that  the  group  velocity  vanishes.  At  the  same  time, 

the  characteristic 

the  group  velocity  is  along 

V  „/V„  ,  =  -2b/(k2dc/dtft)  =  0. 

gk  g  ip 

(69) 

From  this,  we  may  draw  a  physical  picture  of  the  role  of  the  characteristics. 
Plane  waves  emitted  by  a  localized  source,  for  which  k  is  perpendicular  to  a  char- 

acteristic,  do  not  propagate  across  the  characteristic.  Both  the  phase  velocity 
and  group  velocity  vanish,  and  so  the  waves  may  accumulate  along  the  character¬ 
istic.  For  a  sufficiently  steady  situation  (long-lived  source  plus  steady  back¬ 
ground  plasma),  the  waves  may  build  up  to  any  arbitrary  amplitude  within  the  con¬ 
text  of  the  linear  approximation.  Because  the  points  in  space  just  across  the  char¬ 
acteristics  are  not  causally  related  to  those  on  the  side  occupied  by  the  disturb- 
anc  e,  this  disturbance  may  form  a  shock  wave  of  arbitrary  amplitude.  Since  the 
plane  waves  are  supposed  to  extend  an  indefinite  distance  from  the  source,  the 
shock  wave  can  extend  to  large  distances  from  the  source. 


From  the  above  discussion,  we  can  make  clear  what  we  mean  by  "shock-like 
disturbance,  "  a  term  used  throughout  this  paper.  An  observer  passing  through  a 
characteristic  or  resonant  surface  associated  with  a  particular  frequency  would 
measure  a  sudden  jump  in  the  power  at  that  frequency:  the  power  would  be  zero  on 
one  side  and  finite  on  the  other.  However,  since  the  characteristics  themselves 
are  frequency-dependent,  the  disturbance  associated  with  a  band  of  frequencies 
has  a  finite  width:  the  medium  is  dispersive.  Such  a  disturbance  is  quite  unlike 
the  Cherenkov  cone  or  sonic  boom  associated  with  nondispersive  light  waves  or 
sound  waves  (or  fast  mode  waves,  for  that  matter);  these  shocks  show  a  sudden 
jump  in  power  across  a  single  shock  surface  at  all  frequencies. 

Eq.  (69)  will  be  important  in  the  analysis  of  analytic  solutions.  It  can  be 
written  as 

V1W'<*o,«p>‘.,k*  <70> 

where 

*U0€P  *  fp2  '  fH2>/<€0  +  (71) 

From  Eq.  (70),  it  follows  that  when  <  0,  the  group  velocity  is  inside  the  char¬ 
acteristics,  while  when  £  ,  >  0,  the  group  velocity  is  outside  the  characteristics. 
"Inside"  and  "outside"  here  mean  "toward  the  field  direction"  and  "away  from  the 
field  direction,  "  respectively.  This  distinction  will  lead  to  two  separate  classes 
of  shock-like  disturbances  caused  by  spacecraft  corresponding  to  spacecraft  ve¬ 
locity  essentially  across  the  field  and  along  the  field,  respectively.  There  will 
also  be  two  different  ways  for  waves  excited  by  spacecraft  to  spread:  along  the 
field  and  across  the  field,  respectively. 

The  results  of  this  section  can  be  illustrated  in  the  limit  of  low-frequency 
waves  for  which  solutions  can  be  obtained  analytically.  We  state  the  following 
without  proof.  In  the  limit  of  low  frequency,  the  characteristics  are  determined 
by  the  equation 

[d  2  -  <U)2  U>  2>d  2  -  tu2  'V,2J<jB  -0  (72) 

7  gm  x  A  J  y 

where  the  geometric  mean  frequency  (see  Stix,  *  p.  32)  is  defined  by 


The  phase  velocity  is 


ot2/k2  =  VA2  [cosV  -  (t*/2/cjgm2>  sin2#].  (74) 

The  directions  of  the  group  and  phase  velocities  are  related  by 

tan  #g  =  -<a>2/(jgm2)  tan#.  (75) 


The  square  of  the  group  velocity  is 

V  2  5  V,2  [cos2#  -  (tit  2 /CJ2 )  sin2#  ]  /cos4#  .  (76) 

g  A  >-  ^g  gm  g  g 

The  phase  velocity  goes  to  zero  at 

tan#  =  t-  cj  lu)  (77) 

gm 

and  the  group  velocity  goes  to  zero  at 

tan#g  =  (78) 

Iiecause  the  angles  between  the  characteristics  and  field  are 


t  tan  1  (tu/w  ), 
gm 


(79) 


the  main  results  of  this  section  follow  easily.  When  k  is  perpendicular  to  a  char¬ 
acteristic,  the  phase  velocity  is  zero  [Eq.  (74)],  the  group  velocity  is  along  the 
characteristic  [Eq.  (78)]  ,  and  the  group  velocity  is  zero  [Eq.  (76)]  .  The  phase 
and  group  velocities  are  shown  in  Figure  2,  along  with  the  characteristics.  We 
show  (j/Wgm  =  0.2  for  purposes  of  clarity.  In  fact,  the  low-frequency  approxima¬ 
tion  used  requires  w<<  Q .  or  at/iugm  <<  ( mg/ m .)  -  0.  023.  The  magnetic  field 

in  the  figure  is  vertical,  as  indicated  by  the  arrow. 

2 

We  note  that  Eq.  (76)  is  different  from  the  results  illustrated  by  Musielak. 

Ilis  low-frequency  group  velocity  loci  do  not  go  to  zero  at  any  angle. 


I'ii’urc  2.  The  Phase  :in<l  Croup  Velocities  of  Low  Frequency  Plasma 
Waves,  (liven  as  Functions  of  Ancle  With  Respect  to  the  Magnetic 
Field.  The  direction  of  the  magnetic  field  is  indicated  by  the  arrow. 
Ttie  characteristics  are  also  plotted.  The  ancle  between  the  charac¬ 
teristics  and  the  magnetic  field  is  exaggerated  for  clarity. 


I.  I'HOCHFSSIM.  W  \\  K  SOI.I  I  IONS 

With  the  discussion  of  the  characteristics  and  the  propagation  of  discontinui 
ties  as  background,  we  now  turn  to  solutions  of  the  equations.  It  is  well  known 
that  equations  such  as  Kqs.  (25)  and  (26)  or  (36)  can  be  solved  by  the  method  of 
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progressing  waves  (see  Courant  and  Hilbert,'  Chap.  \  I).  For  the  present  pur¬ 
poses,  we  will  find  approximate  solutions  of  Kqs.  (27)  and  (28).  To  do  this,  we 
look  for  solutions  of  the  form 


m. 


[d(,l)(U)a  (fJ.tf)  ■  d{n)(V)<  (fJ,  I/)] 


0  =  L"=  o  [<5(n)<0)bn(Af,  v)  +  6{n\v) dn(V.v)] 


where  U  and  V  are  the  characteristics 

M  =  z  +  ^(«0 1  €p^x  <82) 

v  =  z  -  V u Q/  «p)x  (83) 

[see  Eq.  (41)].  In  Eqs.  (80)  and  (81),  ,5  ^(x)  is  the  Dirac  delta  function,  6^n\x) 
is  its  nth  derivative  when  n  <  0,  and 

<*(n)(x)  =  xn/n ! ,  x  >  0 

=0.  x  <  0  (84) 

when  n  >  0.  Because  of  the  Fourier  transform  in  time,  the  delta  function  indi¬ 
cates  a  source  oscillating  at  frequency  u.  Because  of  the  form  of  the  expansions 
(80)  and  (81),  the  disturbance  is  assumed  to  vanish  outside  the  characteristics. 

The  results  to  be  derived  below  hold  only  for  a  source  at  rest  in  the  proper 
frame  of  the  plasma.  This  restriction  arises  because  of  our  treatment  in  config¬ 
uration  space  rather  than  Fourier  space.  While  the  Doppler  shifts  associated 
with  movinp  sources  are  easily  treated  in  Fourier  space,  the  reversion  to  config¬ 
uration  space  represents  an  intractable  analytical  problem.  It  is  possible  to  treat 
it  numerically,  however.  Remarks  about  moving  sources  are  deferred  until  later 
sections. 

The  treatment  that  we  have  outlined  does  not  properly  treat  the  source  cur¬ 
rents  associated  with  the  spacecraft  interaction  that  cause  the  fluctuations.  To  do 
so  here  would  unduly  complicate  the  presentation.  The  equations  with  source 
terms  are  listed  in  the  Appendix,  along  with  a  discussion  of  the  modifications  re¬ 
quired  to  use  them. 

Substitution  of  Eqs.  (80)  and  (81)  into  Eqs.  (27)  and  (28)  followed  hv  the  equat 
ins  to  zero  of  the  coefficient  of  each  (5*n)(/i)  leads  to  the  recursion  relations 
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(1  +  <0/€p)[b.  +  2(dti  -  +  2(dtl  -  d„)  bj.2] 


£  Pbi  -  2  +  £  Hai  -  2  =  °*  1  >  1 


■“B/'pJb^a^o-^H/V1^-  ^)b0  =  0 

'(<H/£P)bi  +dv  ai-  1  '2(tH/£  P)(S|U  '  dV)bi-  1  +  du  Qvai-2 


£Pai  -  2  £Hbi  -  2  +  ^ft  *  ^  bj  _  2  =  0,  i  >  1 


The  recursion  relations  for  and  d  are  entirely  similar  and  can  be  obtained  by 
interchanging  a  and  c,  b  and  d,  and  fj  and  V.  To  specify  the  solution,  we  choose 


aQ(0,  V )  =  1 


a^O,  V)  =  0.  i  >  0, 


This  choice  implied  that  a  is  a  delta  function  at  fi  =  0,  V  =  0,  The  meaning  of  the 
solution  with  these  initial  conditions  is  that  it  represents  the  manner  in  which  the 
disturbance  that  results  from  an  oscillating  point  source  propagates.  The  solu¬ 
tions  for  these  initial  conditions  are 


ao  =  1 


a  j 


a  g  =  (1/2  )(<„(/)* 


(l/3!)(  -  [(fH2ef,)/(£0  *  (p)]  V 


ho  =  0 


h  =  0 


h2  =  '£1’*U  <£0  ‘  £f' 


Nj-vi' 


r-.  •.  r, < 
*V'  *  ,  * 

vv;«\y; 

»  >  v 


y-y'vy 


■  s  *•  /- 


-  >  .  VJ 

_-*y 
•r. '  * 


■VET 


,V  .*•  n“’  . 


b3  =  ^£P£H^£0  +  £P^  *£**^ 


b4  =  -€H^l/2Hf*U)2  +  (2e0eP  +  3(P2  -  2fH2)/(£0  +  £pd 


£*~(£0£P  +  £P  *€H^(eO+€P)' 


(99) 

HRL  •*< 

(100) 

"*n" 

—  •  m  ^ 

(SW&fr, 


From  these  terms,  it  is  possible  to  describe  the  nature  of  the  solution.  The  delta1 
function 

«=  (I  +  (p/f0)  6 By  -  (fH/e0)  dBx  =  Hft)  (102) 

lies  along  the  characteristic  line  //  =  0,  that  is,  along  z  =  -  p)x.  There  is 

no  corresponding  delta  function  in  /?  -  (1  +  €p/<Q)dBx  +  ^£H^£0^^y*  There  are, 
however,  delta  functions  in  both  6B  and  <5B  . 

y  x 

The  coefficient  a  j  represents  the  jump  in  Of  across  jf  =  0.  This  jump,  which 
is  initially  zero  (because  of  our  choice  of  initial  conditions),  prows  linearly  with 
distance  away  from  the  source.  Again,  there  is  no  jump  in  ft,  but  there  are 
jumps  in  both  and  5P>x.  The  coefficient  a2  represents  the  jump  in  the  first 
derivative  of  a;  it  grows  quadraticall.y  away  from  the  source.  The  jump  in  the 
first  derivative  of  ft  does  not  prow  quadractically  near  the  spacecraft,  however. 
Eq.  (98)  shows  that  the  jump  in  the  first  derivative  of  ft  is  discontinuous,  even  at 
the  origin. 

The  nature  of  the  solutions  is  that  of  coupled  waves.  The  coupling,  as  men¬ 
tioned  in  connection  with  Eqs.  (34)  and  (35),  is  determined  by  €p.  When  tp  Poes 

to  zero,  the  solutions  for  the  a  and  b  are 

n  n 


l-v-v: 


an  =  (€pP)  / n! 


h  =  0. 

n 


Summing  the  series  (80)  and  (81),  we  have 


01  =  2^n"o  «P^>n/<n!)2 


=  2J  [V(-2t  jig.] 


(105) 
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where  Jg(z)  is  the  Bessel  function  and 


£=  0. 


(106) 


The  factor  of  2  in  Eq.  (105)  appears  because  of  the  contributions  of  the  c  in 
Eq.  (80).  Eq.  (105)  is  the  well  known  fundamental  solution  of  Eq.  (34). 

From  this  discussion,  it  appears  that  the  first  few  terms  of  or  are  the  same  as 


a=  J0|V(-2€*^>].  <  0 


(107) 


and  that  modifications  appear  because  of  mode  coupling  in  a^,  bg»  and  higher  order 
terms.  The  former  result  is  easy  to  prove.  By  eliminating  the  low  order  terms 
in  the  recursion  relations,  so  that  they  reduce  to 


■■a 


0."vV\ 

V. 


-  ..'v  sai 


IvVV: 


(1  +f0/eP)bi  +  Vi-2  +  '  •  •  =  0 


a.  ^  j  —  ^H^P^i  *  cpai  -  2  ^  *  *  *  ” 


(108) 

(109) 


one  immediately  obtains 


«*ai-l  =  (*al-2  + 


(110) 


so  that  a  is 


a  -  JQ [V(-2 ^ i/)]  +  ...,«*<  0 


(111) 


where  the  ellipsis  represents  the  remaining  terms  that  arise  because  of  the  cou¬ 
pling. 

In  contrast  to  the  low-frequency  Alfven  waves,  the  solution  for  higher  fre¬ 
quency  waves  depends  on  the  sign  of  £+.  If  >0,  it  would  appear  that 


a=  l0[V(2«%/iv)]+  .  .  .  .  e*  >  0 


(112) 
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where  IQ(z)  is  the  modified  Bessel  function.  The  asymptotic  form  of  Eq.  (112)  for 
laree  argument  is 
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Ct  =  e^ /✓(27 r£)  *  .  .  .  .  c*  >0  (113) 

where  sl  =  -2 1^.(1  V.  which  would  suggest  an  exponential  increase  in  the  wave  ampli¬ 
tude  with  z. 

Such  exponential  behavior  presumably  does  not  occur  for  the  following  reason: 
In  l'.q.  (70),  we  saw  t hat  the  sign  of  determined  the  side  of  the  characteristics 
upon  which  the  waves  could  propagate.  If  e  >  0,  the  waves  propagate  outside 
(away  from  the  field)  of  the  characteristics,  while  the  form  of  Eqs,  (80)  and  (81) 
assumes  that  the  disturbance  propagates  inside  and  vanishes  outside.  The  intui¬ 
tively  correct  result  is  obtained  hy  replacing  the  delta  functions  of  Eqs.  (80)  and 
(81)  by  and  6^n\- V)  if  >  0.  With  an  assumed  solution  of  this  form, 

the  sign  of  //or  V  (but  not  both)  is  flipped  in  each  expansion  with  coefficients  an 

and  c  .  The  result  is  that 
n 

a  --  J  [✓(-2  e.fiv)]  0  (114) 

that  is,  the  solution  is  of  the  same  form  as  Eq.  (Ill)  hut  vanishes  inside  the  char- 
aeteristies. 

The  behavior  of  t in  regions  where  real  characteristics  exist,  is  shown  in 
Figure  3  for  t he*  same  low  Earth  orbit  conditions  as  Figure  1.  Note  that  e*  <  0 
in  the  Alfven  and  whistler  regions  and  >  0  in  the  upper  hybrid  region. 

It  should  lie  noted  that  the  correspondence  of  "inside"  and  "outside"  with  the 
sign  of  was  a  purely  local  result:  it  did  not  determine  whether  the  group  veloc¬ 
ity  was  always  inside  or  outside  of  the  characteristics.  Because  of  the  compli¬ 
cated  nature  of  the  group  velocity  loci,  and  the  fact  that  the  group  velocity  along  or 
across  the  field  is  not  related  to  €  such  a  determination  would  presumably  have 
to  come  from  a  numerical  search.  So  we  are  faced  with  the  following  alternatives: 
(Cither  our  correspondence  is  correct,  or  the  wave  amplitudes  can  grow  exponen¬ 
tially  with  distance  near  the  spacecraft.  While  the  existence  of  growing  waves 
would  be  surprising,  it  is  not  impossible.  Even  though  the  driving  Iq  term  grows 
rapidly  near  the  spacecraft,  the  remaining  terms  could  easily  conspire  to  reduce 
it  at  great  distances.  Because  of  the  hyperbolic  nature  of  the  problem,  there  is  no 
boundary  condition  far  from  the  spacecraft  which  can  be  used  to  eliminate  this  so¬ 
lution.  The  only  apparent  constraint  is  that  the  asymptotic  wave  amplitudes  are 
nondecreasing.  This  is  essentially  an  energy  conservation  argument:  If  the  asymp¬ 
totic  wave  amplitudes  continued  to  grow,  the  steady  source  could  not  provide  suf¬ 
ficient  energy  to  sustain  them.  Nor  could  the  cold  plasma  provide  the  energy,  as 
it  is  known  to  be  stable.  It  would  be  a  worthwhile  exercise  to  look  for  such  waves; 


.  .."  V 


23 


4«, 


iff, 


a>  (SEC) 

Figure  3.  The  Parameter  €*  Plotted  as  a  Function  of  Wave  Frequency.  The 
properties  of  the  plasma  are  the  same  as  Figure  1.  In  the  Alfven  region  (o»< 

2  x  10®  sec  “1)  and  the  whistler  region  (10®  <  <  7  x  10®)  f  *  is  negative.  In  the 

upper  hybrid  region  (8  x  10®  <  U>)  t  if  is  positive 


they  could  provide  a  source  of  extremely  rapid  spatial  wave  growth  near  the  space¬ 
craft. 

Several  remarks  should  be  made  about  the  nature  of  these  solutions. 

(1)  Despite  the  fact  that  the  solutions  are  analytic,  thev  can  only  be  regarded 
as  qualitative;  they  are  presented  only  for  the  purpose  of  identifying  typical  physi¬ 
cal  features. 

(2)  The  solutions  are  only  two-dimensional;  solutions  can  change  significantly 
when  the  dimensionality  of  the  problem  is  chaneed.  The  best  known  example  is 


v'.’S'.' 


light:  three-dimensional  solutions  are  delta  functions,  while  two-dimensional  solu¬ 
tions  are  Bessel  functions.  While  both  have  sharp  wave  fronts,  the  two-dimen¬ 
sional  solutions  have  trailing  wakes  entirely  missing  in  three  dimensions. 

(3)  The  solutions  of  this  section  are  for  delta  function  disturbances  in  the 
magnetic  field.  A  better  physical  description  would  be  the  solutions  for  a  delta 
function  current  source.  (For  further  details,  see  the  Appendix. ) 

(4)  The  squares  of  the  results  of  this  section  represent  the  spectral  power 
density  of  the  fluctuations  resulting  from  a  two-dimensional  source  stationary  in 
the  proper  frame.  For  applications,  the  power  spectral  density  should  be  calcu¬ 
lated  for  specified  current  sources  (see  Appendix). 


5.  STEADY  PLASMA  STR1T.T1  RK  AROI'ND  A  MOVING  SOI  RCK 

The  present  section  is  somewhat  outside  the  main  line  of  development  and  can 
be  omitted  on  a  first  reading.  The  purpose  is  to  develop  a  formal  method  for  cal¬ 
culating  the  steady  plasma  structure  around  a  moving  source.  Because  of  the  in¬ 
tractable  analytic  nature  of  the  result,  we  have  computed  results  only  in  two 
simple  cases.  However,  these  cases  are  instructive  and  may  help  the  uninitiated 
reader  to  make  contact  with  the  subject.  The  formal  results  include  the  entire 
steady  radiation  field,  which  has  contributions  from  waves  of  all  frequencies.  As 
such,  it  should  be  able  to  serve  as  a  useful  starting  point  in  future  investigations. 
In  contrast  to  the  rest  of  this  paper,  the  results  of  this  section  are  fully  three-di¬ 
mensional. 

The  problem  that  we  wish  to  investigate  is  the  plasma  response  to  a  source  of 

♦ 

current  moving  with  constant  velocity,  U.  Let  the  source  be  characterized  by  a 
charge  density,  T)  ,  so  that  a  point  source  can  be  represented  by 

J  =  T) U  =  \j<5(x  -  Ut).  (115) 

We  wish  to  use  this  source  current  to  calculate  a  Green's  function,  and  so  have 
chosen  T)  =  1.  The  Ampere-Faraday  Law  is 

Vx  (V  X  <5E)  =  *i;j0E  -  (4TTio//c2)7  (116) 

where  the  represent  the  coefficients  defined  in  Eq.  (10).  Without  working  out 
the  details,  one  can  see  that  the  oeneral  form  of  the  spatial  Fourier  transform  of 
Eq.  (116)  is 

A  . t  (5E  =  (47Tiw/c2)J. 


(117) 
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where  A  is  a  matrix  whose  elements  depend  onT  and  f  „(<*>).  The  solution  for 


dEi  =  (47Titj/c2)Ai."1J 


where  the  inverse  matrix  A  is  given  by 


A  . j '  =  ay(w.  k)/D(<i>,  k). 


Here,  the  a.^(o»,  k)  are  coefficients  resulting  from  the  inversion;  they  depend  on  u 
throuch  the  <^(ot).  D  (<j,  k)  is  the  dispersion  function;  the  dispersion  relation  is 
D(ut,  k)  =  0.  The  Fourier  transform  of  the  Green's  function  for  this  point  source 


G.(w,  k)  =  (47Tiw/c^)[aij«*>,  k)/D(u>,  k)]j.(tj,  k) 


where  JAW,  k)  is  the  Fourier  transform  of  the  current  density  given  by  Eq.  (115) 
and  is  given  explicitly  by 


J.(w,  k)  =  [l./<27r)3]<5(u  -  k.U). 


Then  the  Green's  function  is 


G^t.  t)  =  /d3kd<XT(W.  k). 


After  some  easy  steps,  the  Green's  function  can  be  written  symbolically  in  the 


GAx  -  Tl)  =  0. /d'\  [exp(ik-  U  -  Ft))]  /D(k-f,  IT)  (123) 

where  0.  is  an  operator  defined  bv 

0.  =  [ri/(2/r2c2)](lT*?)ai.(-ir.^,  -i7>.  (124) 

Note  t  t.it  a .  (  depends  on  T*7  through  the  dependence  of  the  onu>.  It  is  a  compli- 
c.'i.  d  function  of  differential  operators.  Its  evaluation  normally  involves  power 
series  of  operators. 

Several  remarks  about  Eq.  (124)  are  appropriate: 
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(1)  The  wave  pattern  is  stationary  in  the  frame  in  which  the  source  is  at  rest. 

(2)  There  are  two  principal  pieces  to  the  Green's  function:  the  operator  0.  and 
the  radiation  integral 

R  —  /dgk[exp(ik*(x  -  Lit))] /D(k«  U.  k).  (125) 

(3)  The  radiation  integral  carries  the  fundamental  information  about  the  ex¬ 
pansion  of  the  waver  phase  speeds,  group  speeds,  characteristics,  and  disturbance 
topology. 

(4)  The  details  of  each  wave  amplitude  are  determined  by  the  operator  0,. 

(5)  R  is  simply  the  spatial  Fourier  transform  of  the  inverse  dispersion  func¬ 
tion. 

(6)  R  uepends  on  all  wave  frequencies  because  of  the  Doppler  shift,  repre- 
sented  through  the  dependence  of  D  on  k*U. 

Eq.  (123)  is  difficult  or  impossible  to  evaluate  in  any  general  way.  We  would 
like  to  illustrate  its  use  through  a  simple  evaluation  of  R  for  two  special  cases: 

(1)  fast  mode  waves,  and  (2)  MHD  Alfven  waves.  To  do  so,  we  shall  drop  all  con¬ 
stant  factors  and  derivatives,  which  are  inessential. 

(1)  The  fast  mode  dispersion  relation  is 

D((j.  kj  =  (J2  -  k2VA2  (126) 

so  that 

D(k.U.  "kj  =  (k.U)2  -  k2VA2.  (127) 

To  proceed,  let  L  =  (Ux.  0,  0)  so  that  the  velocity  is  perpendicular  to  the  field. 
Then 

D(k.U.  kj  -  k  2  -  (k  2  +  k  2)/(U2/V  2  -  1).  (128) 

x  y  Z  n 

in  carrying  out  the  transforms  in  Eq.  (123),  there  are  two  cases  to  consider.  If 
2  2 

U  /V.  >  1,  D  has  two  real  roots  for  k  .  Carrying  out  the  transform  in  the 

A  X 

usual  way,  we  obtain 

R  ~  l//[x2/(U2/VA2  -  1)  -  y2  -  z2]  if  y2  f  z2  <  x2/(U2/VA2  -  1) 


=0,  otherwise 


This  represents  a  fast  mode  shock  for  the  case  that  the  source  speed  exceeds  the 

Alfven  speed.  The  disturbance  vanishes  outside  the  shock  cone,  shows  a  sudden 

jump  at  the  shock  cone,  and  is  nonzero  inside.  This  is  just  Cherenkov  radiation 

in  a  cold  plasma. 

2 

If  L:  /V.  “  <  1,  there  are  two  imaginary  roots  of  k  .  In  this  case,  the  radia- 

A  X 

tion  integral  is 

R  ~  1/  [x2/(l  -  1 :2/Va2)  *  y2  *  z2].  (130) 

This  has  the  same  form  as  Eq.  (129)  except  that  there  is  no  shock  cone. 

(2)  The  MUD  Alfven  mode  dispersion  relation  is 

D(W.  T)  =  U)2  -  (k.  V.)2  (131) 

so  that 

IXT.T.T)  <k«TT)2  -  <k.TA)2.  (132) 

With  the  velocity  again  in  the  x -direction 

OUT.  tr,  T)  kx2l'2  -  kz2vA2.  (133) 

Evaluation  of  the  integrals  then  leads  to 

R  ~  d(y)[d(x  >  Uz./VA)  -  <Mx  -  Uz/VA)].  (134) 

This  solution  is  just  two  delta  functions  in  opposite  directions  along  the  field:  they 

are  the  Alfven  wings  for  a  point  source.  They  are  tilted  by  an  angle  tan  *(U/V.) 

3  A 

with  respect  to  the  field  and  lag  behind  the  source. 

(,.  \ITI.I<  Vl  lovs  TO  SI*  \CK(  H  \FT 

1'p  to  this  point,  we  have  developed  several  apparently  disparate  approaches 
to  our  problem.  It  is  time  to  try  to  unify  them  into  one  coherent  picture.  The 
main  difficulty  we  face  is  to  incorporate  the  motion  of  the  spacecraft  into  our  re¬ 
sults.  Aside  from  the  last  section,  the  results  are  valid  only  for  sources  which 


are  stationary  in  the  proper  frame  of  the  plasma.  For  such  a  motionless  source, 
the  methods  outlined  above  should  provide  a  detailed  accounting  of  the  behavior  of 
waves  and  disturbances  propagating  away  from  the  source. 

What  is  required  is  a  generalization  of  the  characteristics  that  applies  to  the 
wave  field  of  moving  sources.  These  characteristics  should  be  the  carriers  of 
discontinuities  that  originate  at  the  spacecraft.  The  spatial  methods  that  have 
been  developed  above  are  not  particularly  helpful  for  the  following  reasons:  The 
characteristics  in  the  proper  frame  do  not  depend  on  the  group  velocity  because, 
as  we  have  seen,  it  becomes  zero  at  the  characteristics.  On  the  other  hand,  the 
spatial  characteristics  in  a  flowing  plasma  must  have  some  relation  to  the  group 
velocity.  Another  way  of  looking  at  it  is  that  the  angle  of  the  characteristics  in 
the  frame  of  the  moving  source  must  be  determined  by  the  ratio  of  some  speed  in 
the  problem  to  the  source  speed,  U.  This  is  true  in  any  frame;  if  the  source  is 
stationary  in  the  proper  frame,  the  source  velocity  and  the  group  velocity  associ¬ 
ated  with  the  characteristics  are  both  zero. 

The  appropriate  generalization  has  been  known  for  a  long  time  (see  Jeffrey 
6 

and  Taniuti,  p.  186).  If  V  -  U  is  tangent  to  the  group  velocity  surface,  then  the 
♦  «.  S 

direction  -  U  is  characteristic: 

(Vg-U)-F0=O  (135) 

that  is.  the  convected  group  velocity  is  parallel  to  the  characteristic  surface.  The 
construction  of  the  characteristics  is  shown  in  Figure  4;  we  have  used  the  group 
velocity  of  Figure  2  in  this  example.  The  magnitude  of  the  spacecraft  velocity  has 
been  greatly  exaggerated  for  clarity. 

The  characteristics  determined  in  this  way  satisfy  the  usual  requirement  that 
characteristic  surfaces  be  the  carriers  of  discontinuities.  This  is  easy  to  see 
from  the  significance  of  the  group  velocity.  The  characteristics  here  are  just  the 
envelope  of  the  group  velocity  and  so  determine  the  limits  past  which  a  localized 
fluctuation  cannot  spread:  inside  the  characteristic  surface,  the  disturbance  does 
not  vanish,  while  outside,  it  does. 

g 

For  further  details,  see  Jeffrey  and  Taniuti.  For  the  present  purposes,  it 
is  sufficient  to  note  the  following:  Spacecraft  velocities  are  about  8  km/ sec  or  less 
while  typical  group  velocities  are  on  the  order  of  the  Alfven  speed,  several  hun¬ 
dred  kilometers  per  second  or  more.  Unless  there  are  group  velocities  less  than 
the  spacecraft  velocity,  then  no  new  qualitative  features  (that  is,  new  characteris¬ 
tic  surfaces)  will  occur.  In  that  case,  the  analysis  of  this  paper  can  be  used  to 

6.  Jeffrey,  A.,  and  Taniuiti,  T.  (1964)  Non-Linear  Wave  Propagation,  Academic 
Press,  New  York. 


Figure  4.  Construction  of  the  Characteristics  for  a  Moving  Source.  The 
polar  plot  is  the  group  velocity  locus;  the  magnetic  field  direction  is  indi¬ 
cated  by  the  arrow.  The  source  velocity  is  indicated  by  U 

determine  all  of  the  shock-like  disturbances.  It  should  be  emphasized,  however, 
that  larger  source  velocities  would  introduce  qualitatively  new  characteristic  sur 
faces  and  the  associated  shock-like  disturbances;  an  example  is  the  fast  mode 
Cherenkov  cone. 

In  the  limit  of  small  spacecraft  velocity,  it  is  easy  to  see  the  qualitative  na¬ 
ture  of  the  changes  in  the  characteristic  surfaces  and  to  discuss  the  nature  of  the 


spacecraft  interaction  associated  with  them.  Let  us  discuss  the  qualitative  nature 
of  the  changes  for  two  cases:  £*  <  0  (Case  1)  and  e*  >  0  (Case  2). 

Case  1:  If  the  spacecraft  velocity  is  across  the  field,  shock-like  disturbances 
can  occur.  Ahead  of  the  spacecraft,  the  changes  will  be  small.  There  will  still 
be  a  well  defined  characteristic  surface  past  which  disturbances  of  a  given  fre¬ 
quency  cannot  spread.  The  angle  this  surface  makes  with  the  magnetic  field  will 
still  be  given  to  high  accuracy  by  Eq.  (44).  Behind  the  spacecraft,  the  changes 
are  more  significant.  Because  the  waves  with  zero  group  velocity  in  the  proper 
frame  do  not  propagate  through  the  plasma,  once  they  are  produced,  they  will  form 
a  trailing  wake  behind  the  spacecraft.  This  is  simply  a  statement  of  the  fact  that 
there  is  no  characteristic  surface  in  the  spacecraft  wake;  unlike  the  case  of  a  sta¬ 
tionary  source,  there  are  no  solutions  in  the  wake  of  Eq.  (133)  for  n,  as  can  be 
seen  from  Figure  4.  Even  though  there  are  no  trailing  characteristics,  the  solu¬ 
tions  obtained  ( Eqs,  (92)  through  (100)]  should  still  be  reasonably  accurate,  ex¬ 
cept  that  the  trailing  disturbance  will  approach  zero  rapidly  but  not  vanish. 

If  the  spacecraft  velocity  lies  within  the  characteristics,  no  shock-like  dis¬ 
turbances  are  possible.  For  example,  if  the  velocity  is  parallel  to  the  field,  the 
disturbance  propagates  (in  the  spacecraft  frame)  ahead  of  the  spacecraft  at  several 
hundred  kilometers  per  second  and  behind  the  spacecraft  at  the  orbital  velocity 
(several  kilometers  per  second). 

While  the  formal  results  have  been  obtained  in  two  dimensions,  we  may  apply 
them  in  a  qualitative  fashion  to  three  dimensions.  The  application  follows  easily 
from  the  axisymmetric  nature  of  the  actual  waves  around  the  field  in  the  proper 
frame.  There  is  a  two-lobed  disturbance  field.  The  dividing  line  between  the  two 
lobes  trails  the  spacecraft  antiparallel  to  the  spacecraft  velocity.  The  leading 
edge  of  the  disturbance  (ahead  of  the  spacecraft)  is  a  cone  whose  angle  with  re¬ 
spect  to  the  field  is  given  approximately  by  Eq.  (44).  Toward  the  sides  of  the 
cone,  the  disturbance  is  still  sharply  defined,  with  much  the  same  cone  angle. 
However,  the  trailing  disturbance  extends  to  infinity,  growing  narrower  at  large 
distances.  A  qualitative  representation  is  shown  in  Figure  5.  For  the  conditions 
of  Figure  1,  this  configuration  applies  to  Alfven  and  whistler  disturbances. 

Case  2:  If  the  spacecraft  velocity  is  along  the  field,  shock -like  disturbances 
can  occur.  Inside  the  characteristics  ahead  of  the  spacecraft,  there  will  be  no 
disturbance.  The  main  disturbance  itself  consists  of  two  lobes,  associated  with 
the  group  velocity,  directed  across  the  field.  Again,  there  will  be  a  trailing  wake. 
If  the  velocity  lies  outside  the  characteristics,  then  there  are  no  shock-like  dis¬ 
turbances. 

In  three  dimensions,  the  disturbance  is  obtained  bv  rotating  the  two-dimen¬ 
sional  result  around  the  field.  The  result  is  a  conical  leading  disturbance  with  a 
trailing  wave  field.  The  main  disturbance  is  expected  to  lie  between  two  mirror 


symmetry  cones  (for  the  case  that  Li  is  parallel  to  B)  with  a  smaller  amplitude 
disturbance  extending  to  infinity  in  the  wake.  A  qualitative  representation  is 
shown  in  Figure  6.  For  the  conditions  of  Figure  1,  this  configuration  applies  to 
upper  hybrid  disturbances. 


Figure  5.  Qualitative  View  of  the  Three-Dimensional  Limits 
of  the  Wave  Field  for  a  Given  Frequency  and  €  *  <  0.  For 
the  plasma  of  Figure  1,  this  behavior  is  typical  of  the  Alfven 
and  whistler  regions 


Figure  6.  Qualitative  View  of  the  Three-Dimensional  Limits 
of  the  Wave  Field  for  a  Given  Frequency  and  c ^  >  0.  For 
the  plasma  of  Figure  I,  this  behavior  is  typical  of  the  upper 
hybrid  region 

7.  DISCISSION 

We  have  investigated  the  properties  of  waves  that  are  assumed  to  be  produced 
as  the  result  of  the  interaction  of  a  spacecraft  with  a  cold  plasma.  We  have  found 
that  there  are  three  distinct  frequency  regimes  for  which  the  waves  can  propagate 
to  large  distances  from  the  spacecraft  and  cause  shock-like  disturbances.  These 
waves  may  be  roughly  classified  as  Alfven  waves,  whistlers,  and  upper  hybrid 
waves.  In  the  examples  given  in  the  paper,  the  Alfven  and  whist’er  waves  are 
guided  along  the  geomagnetic  field,  while  the  upper  hybrid  waves  propagate  across 
the  field. 

The  question  of  the  production  of  these  waves  has  not  been  dealt  with.  For  the 


purposes  of  this  work,  we  have  simply  assumed  the  existence  of  a  wave  source 
and  have  treated  the  subsequent  evolution  of  the  waves.  The  production  of  Alfven 
wings  is  reasonably  well  understood.  They  result  from  the  steady  charge  separa¬ 
tion  required  to  maintain  zero  electric  field,  as  seen  by  a  comoving  observer,  in¬ 
side  a  conductor  moving  across  a  magnetic  field.  The  resulting  electric  field  at 
the  surface  of  the  conductor  then  propagates  along  the  Alfven  characteristics, 
forming  steady  Alfven  wings. 

The  production  of  higher  frequency  waves  is  not  so  easy  to  understand  in  any 
detail.  In  general,  the  wave  spectrum  associated  with  the  spacecraft  interaction 
should  extend  up  to  a  frequency 

ot^U/L  (136) 

where  L  is  a  typical  dimension  of  the  spacecraft.  For  a  spacecraft  about  10  m 

5  - 1 

across,  the  maximum  frequency  is  about  10  sec  .  If  smaller  scale  structures, 
such  as  wing  tips  and  bay  door  edges  (or,  in  general  any  structure  with  a  small 
radius  of  curvature),  are  considered,  then  the  maximum  frequency  is  raised  per¬ 
haps  two  orders  of  magnitude.  So  it  appears  that  wave  frequencies  as  high  as  the 
whistler  and  upper  hybrid  regions  are  possible. 

On  the  other  hand,  the  formal  result  for  cold  plasmas  is  that  the  spacecraft 
interaction  produces  only  zero  frequency  waves  as  observed  in  the  frame  in  which 
the  spacecraft  is  at  rest,  as  mentioned  in  the  discussion  following  Eq.  (123). 

While  there  are  contributions  to  the  steady  structure  from  waves  of  all  frequencies 
[see  Eq.  (123)],  it  is  not  at  all  clear  that  waves  will  be  produced  that  oscillate  in 
the  spacecraft  frame. 

This  makes  clear  an  underlying  assumption  of  the  paper:  While  the  propagation 
of  cold  plasma  waves  has  been  treated  in  some  detail,  we  rely  on  processes  out¬ 
side  of  cold  plasma  physics  to  actually  produce  the  waves.  While  little  is  known 
theoretically  about  such  wave  production,  there  are  two  lines  of  attack.  One  is  to 
measure  the  wave  field  near  the  spacecraft.  Given  a  measured  source,  the  pre¬ 
sent  formalism  describes  the  subsequent  evolution  of  the  waves.  The  other  is  to 
generate  a  self-consistent  model  of  the  steady  state  Shuttle  environment  and  ana¬ 
lyze  it  for  dominant  plasma  instabilities  and  again  compute  the  evolution.  Both 
types  of  effort  are  in  early  stages  of  development. 


34 


'  J  'J  -T  * 


•XtV'. 


2. 

3. 

4. 

5. 

6. 


.  n.  me  i  neory  oi  riaarna  waves,  ivicuraw-miii  new  i 

Musielak,  Z.E.  (1984)  M.I.  T.  Center  for  Space  Research  Report  No. 

CSR-TR  -84-3. 

Drell,  S.D. .  Foley,  H.  M. ,  and  Ruderman,  M.A.  (1965)  Drag  and  propulsion 
of  large  satellites  in  the  ionosphere:  An  Alfven  propulsion  engine  in  space, 

J.  Geophys,  Res.  70:3131. 

Barnet,  A. ,  and  Olbert,  S.  (1985)  Radiation  of  plasma  waves  by  a  conducting 
body  moving  through  a  magnetized  plasma,  submitted  to  J.  Geophys.  Res. 

Courant,  R.,  and  Hilbert,  D.  (1962)  Methods  of  Mathematical  Physics  (Vol.  II) 
Interscience.  New  York. 

Jeffrey,  A.,  and  Taniuiti,  T.  (1964)  Non-Linear  Wave  Propagation,  Academic 
Press,  New  York. 


Appendix 


The  source  terms  for  the  wave  equations  can  be  introduced  in  the  usual  way. 
The  Ampere -Faraday  Law  is 

x  dE)  =  +  (47ri<j/c2)J  (Al) 

where  J  is  the  source  current.  The  equations  for  the  field  fluctuations  are 

[az2  +  <«p/« 0>dx2  -  ep]^By  «  ‘h4Bx  '  (47r/c)[^zJx  +  <fP/t0)dxJz)  (A2> 

<ex2  +  Qz  *  €P)4Bx  =  -(H(1  +  6o’ldx2)4By  +  [4n,c)dzJy  (A3) 

The  equations  for  the  amplitudes  a  and  /3  are 

[az  ‘  (eo^£P)0x  "  €p]ct  =  CH(1  *  f0  &x 

+  (47T/C)  [-  (1  -€0/ep)azJx  -  (€H/€0)5zJy  -  (1  *«p/«o)(ep/€0>5XJZ^  (A4) 

<az2  -  ax2  -  fp)^=  -«Ha 

M47T/c)[-«H/€0)azJx  *  (1  -  VWy  •  (ePeH/e02)axJz^I  <A5) 


In  the  text,  the  right-hand  sides  were  treated  as  delta  functions.  A  proper  physi¬ 
cal  treatment  requires  that  the  current  density  be  the  delta  function  source.  The 
qualitative  result  is  that  the  expressions  in  the  text  must  be  differentiated  with  re¬ 
spect  to  the  appropriate  variable  (for  example,  z  if  the  source  term  is  d  J  )  and 
then  multiplied  bv  the  appropriate  coefficient  [for  example,  -(47T/c)(l/€^/fp)  for 
the  leading  term  in  Eq.  <A4 )]. 


